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1.  INTRODUCTION 


In  this  first  of  two  papers  we  report  the  a  priori  calculation  of  the  electric  dipole  transition 
moment  as  a  function  of  the  C-C  bondlength  for  the  C^X  electronic  transition  in  acetylene.  This 
moment  is  then  used  to  predict  the  oscillator  strength  for  absorption  and  the  electronic  magnetic 
circular  dichroism  (MCD)  spectra  using  a  newly  implemented  code.  The  0“X  electronic  transition 
occurs  in  the  153-140  nm  region,  with  the  VjCO-O)  (C-C  synunetric  stretch)  at  151.8  nm  (65860  cm  *) 
(Price  1935;  Wilkinson  1958;  Nakayama  and  Watanabe  1964;  Moe  and  E>uncan  1952;  Gadanken  and 
Schnepp  1976).  Previous  experiments  (Price  1935;  Wilkinson  1958;  Nakayama  and  Watanabe  1964; 
Moe  and  Duncan  1952),  as  well  as  ab  initio  calculations  (Peric,  Buenker,  and  Peyerimhoff  1984; 
Kanuner  1974;  Demoulin  and  Jungen  1974;  Demoulin  1975),  have  assigned  the  C  state  to  the  linear 
Rydberg  l'n„  state. 

Demoulin  and  Jungen  (1974)  used  an  ab  initio  quantum  chemical  approach,  but  only  at  the  SCF 
level,  to  predict  an  absorption  oscillator  strength  of  0.072  for  the  (vertical)  0“X  transition.  An  early 
experimental  estimate  of  the  oscillator  strength  was  reported  by  Moe  and  Duncan  (1952)  to  be  0.062, 
while  a  mote  recent  study  by  Gedanken  and  Schnepp  (1976)  reports  both  the  MCD  and  absorption 
spectra.  Unfortunately,  in  the  absorption  spectrum,  the  well  separated  band  maxima  appear  to  reside 
upon  a  very  broad  band  of  non-negligible  intensity.  No  attempt  was  made  to  reveal  the  source  of  this 
broad  band;  however,  an  attempt  was  made  to  subtract  out  the  intensity  of  this  underlying  band  from 
that  which  was  assumed  to  belong  only  to  the  0~X  transition.  After  having  done  this,  they  report  the 
values  for  the  square  of  the  electric  transition  dipole  moments  for  the  first  four  vibrational  bands. 

From  this,  one  can  readily  obtain  the  oscillator  strengths  for  absorption  which  sum  to  0.110  for  these 
four  vibrational  bands. 

The  MCD  spectrum  shows  a  well  resolved  vibrational  structure  corresponding  to  Vj,  with  a 
dispersive  band  shape  indicative  of  a  degenerate  C  state.  In  contrast  to  the  absorption  spectrum,  the 
MCD  shows  no  obvious  signs  of  overlapping  with  the  broad  band  seen  in  absorption.  If  one  then 
extracts  the  electric  transition  dipole  moment  from  the  MCD  intensities  (based  on  a  magnetic  dipole 
moment  of  0.5),  one  gets  the  value  of  0.11  for  the  sum  over  the  four  vibrational  bands,  which  agrees 
well  with  the  value  extracted  from  their  "corrected"  absorption  intensities.  Thus,  the  estimates  of  the 
oscillator  strength  vary  over  a  wide  range  from  0.062  to  0.11,  suggesting  the  need  for  an  estimate  of 
this  transition  intensity  based  upon  high  quality  ab  initio  calculations. 
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A  transition  from  a  non-degenerate  ground  state  to  a  degenerate  excited  state,  as  in  the  case  in  the 
C^X  transition,  should  have  a  large  contribution  to  the  MCD  from  the  ;?-term  (Stephens  1976;  Piepho 
and  Schatz  1985).  In  this  study,  it  is  assumed  that  the  entire  MCD  intensity  can  be  described  by  this 
J5-term,  which  is  calculated  via  ab  initio  quantum  chemical  methods.  The  wavefunctions  will  be 
generated  using  a  large  atomic  gaussian  basis  set  and  state-averaged  MCSCF  plus  configuration 
interaction  (Cl).  In  the  second  paper,  we  will  extend  the  calculations  (both  oscillator  strengths  and 
MCD  intensities)  to  the  190-170  nm  spectral  region  which  includes  the  B<~X  transition.  An  excellent 
review  of  the  literature  pertaining  to  the  C*-X  transition  (as  well  as  other  singlet  states)  can  be  found 
in  the  paper  by  Peric,  Buenker,  and  Peyerimhoff  (1984). 

2.  METHODS 

2.1  Quantum  Chemical  Approach.  The  Cl  method  used  to  obtain  the  zeroth-order  wavefunctions 
is  the  symbolic  matrix  element,  direct-CI  method  of  Liu  and  Yoshimine  (1981).  The  molecular 
orbitals  (MO)  used  as  a  basis  set  for  the  Cl  expansions  were  obtained  from  a  state-averaged,  multi- 
configuration  self-consistent  field  complete  active  space  (SA-MCSCF)  approach  (Diffenderfer  and 
Yaricony  1982).  The  nonactive  MOs  from  the  SA-CASSCF  were  rotated  to  be  eigenfunctions  of  the 
closed  shell  Fock  operator  corresponding  to  the  state-averaged  density.  The  active  molecular  orbitals 
are  taken  as  the  natural  orbitals  determined  from  the  MCSCF.  The  SA-CASSCF  procedure  is  the 
general  second-order,  density  matrix-driven  MCSCF  algorithm  of  Lengsfield  (1980).  From  the 
SA-CASSCF  procedure,  one  obtains  a  set  of  molecular  orbitals  determined  by  minimizing  the  energy 
functional 

E  (1) 

where  the  ’s  are  the  eigenfunctions  of  the  H®,  the  non-relativistic  Hamiltonian  operator,  in  the  space 
of  the  CASSCF  expansion 

'F:-  Ec/vi;,.  (2) 

I 

Above,  the  are  the  nonnegative  weighting  factors  for  the  electronic  states  which  do  not  vary  as  a 
function  of  intemuclear  separation,  and  the  configuration  state  functions  (CSF)  composed  of  the 
MOs  determined  in  the  SA-CASSCF.  The  weights  were  chosen  to  provide  a  balanced  description  of 
the  states  of  interest. 
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2.2  Spectral  Intensities.  The  oscillator  strength  for  the  electronic  transition  from  the  ground  state 
X(l'Zp  to  the  excited  state  CCl’llu)  is  given  by 

/=lA£|<X|/?r"|C>P  (3) 

with  the  transition  energy,  AE,  given  in  au’s.  In  the  reduced  spacial  symmetry  of  this  becomes 


/=  i  A  £  (  I  <X(A,)  1  ^1 C  ( r^)>  r  +  I  <X  (A,)  I  ^1 C  (n,^)>  p  ) .  (4) 


The  MCD  intensity  for  a  transition  involving  a  non-degenerate  ground  state,  X,  and  a  possibly 
degenerate  excited  state,  I ,  is  given  by  the  difference  between  the  decadic  molar  extinction  coefficients 
for  left  and  right  circularly  polarized  light,  i.e.,  Ae  =  Cl  -Er  (Stephens  1976;  Piepho  and  Schatz  1983) 


Ae 

T 


653.24, 


fa/ 


•H 


(5) 


where  E  is  energy,  f  a  lineshape  function,  H  the  magnetic  field,  and  J?,  is  given  by 
a;  -  ^  1  I  ?|  .  >•  (<X  I  *  <4 . 1  X> )  ] 


(6) 


with  4  and  m  being  the  magnetic  and  electric  dipole  moment  operators,  4g  the  Bohr  magneton 
(.4669xlO'^cm'VGauss),  and  X  specifies  the  degenerate  components  of  I.  Equation  6  takes  into 
account  the  fact  that  the  molecules  are  randomly  oriented.  The  units  are  cgs  with  energies  in  cm■^ 
electric  dipole  moments  in  Debye,  and  the  magnetic  field  in  Gauss.  The  operator,  im  is  chosen  in  the 
length  form.  The  symbols  for  the  electric  and  magnetic  dipole  moment  operators  were  chosen  to  be 
consistent  with  the  literature  on  MCD  published  by  Stephens  (1976)  and  co-woikers. 

From  Equation  6,  one  sees  that  the  term  is  derived  for  a  truly  degenerate,  excited  state,  I ,  and 
that  a  MCD  band  attributed  to  such  a  state  will  exhibit  a  characteristic  derivative  lineshape.  The  band 
shape  is  often  used  to  determine  if  a  transition  involves  a  degenerate  excited  state.  The  reader  is 
referred  to  Stephens  (1976)  and  Piepho  and  Schatz  (1983)  for  detailed  discussions  of  MCD  theory, 
applications,  and  band  shapes. 
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As  noted  earlier,  the  C^X  band  system  shows  an  essentially  pure  derivative  band  shape 
consistent  with  a  linear,  degenerate  upper  electronic  state.  Based  on  the  assumption  that  the 
intensity  of  the  MCD  band  comes  strictly  from  the  term,  we  get  (in  symmetry) 


21m 

~T 


:C(n„)  |M,iqn,)><X(A,)  |mJC(n„)><qi;,) 


(7) 


with  the  factor  of  2  replacing  the  summation  over  the  degenerate  components.  The  prime  on 
signifies  that  the  matrix  elements  involve  only  the  electronic  part  of  the  wavefimctions. 


In  order  to  relate  the  observed  spectra  to  our  theoretical  results,  we  use  the  method  of  moments  as 
applied  to  MCD  and  described  by  Stephens,  Mowery,  and  Schatz  (1971).  The  n***  moment  of  a  MCD 
band  can  be  defmed  as 

<Ae>,-  f.£i(£-£")"d£.  (8) 

where  E°  is  an  average  energy  defmed  by  <A£>i=0.  For  the  acetylene  transition  being  studied  here, 
the  first  moment  is 

<A  e>]  =  .01525^  (9) 

with  defined  in  Equation  7. 

Finally,  we  need  to  obtain  the  vibrationally-averaged  electric  transition  dipole  moments  between 
the  two  bound  electron  states  in  Equation  7.  For  this  purpose,  a  section  of  the  potential  energy  surface 
(PES)  was  calculated  which  approximates  the  normal  mode  describing  the  C-C  symmetric  stretch,  with 
the  vibrational  coordinate,  "R,"  being  the  C-C  bond  length.  The  C-H  bond  lengths  are  held  fixed  and 
the  molecule  remains  linear.  Using  these  PESs  to  calculate  the  vibrational  wavefimctions  Xv"(^)  ^nd 
Xv'(R)  for  the  'F(X)  and  ^(Q  electronic  wavefimctions,  respectively,  we 

5(;f;C):.v  =  <X.«(/?)|<^(^)|'«JC(nj>|x..(R)>  (10a) 

5(C;X),V  =  <x,.(R)|<C(n,)|mjX(A,)>|X,..(/?)>.  (10b) 

The  XvWs  are  obtained  by  numerically  solving  the  radial  Schroedinger  equation  for  nuclear  motion 
while  ignoring  rotational  effects.  The  integrals  in  Equations  10a  and  10b  are  then  solved  numerically 
using  Simpson’s  rule. 
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(11) 


The  final  forms  for  f  and  Ai,  which  includes  vibrational  averaging,  are 

iA£{  •  is(e;4’vP) 

^.(.■v,-^(<c(q„)lml<:(n„)>-5(x;(:);.v  •«(<?;#):.,..).  (12) 

Note  that  the  magnetic  dipole  transition  moment  in  Equation  12  is  a  constant  (i.e.,  one-half  the  orbital 
angular  momentum  quantum  number  for  a  11  state)  with  respect  to  the  R(C-C)  stretch  and,  therefore, 
excluded  from  the  vibrational  treatment. 

3.  DETAILS  OF  CALCULATIONS 

The  basis  set  used  throughout  this  study  consists  of  the  primitives  (lls.7p)  contracted  to  [5s.4p] 
on  the  carbons  and  (6s)  contracted  to  [3s]  on  the  hydrogens  (McLean  and  Chandler,  unpublished 
results).  The  carbon  atoms  were  augmented  with  two  uncontracted  d-type  polarization  functions  arid 
the  hydrogens  each  received  one  uncontracted  p-type  ftmction.  To  the  center  of  the  C-C  bond  was 
added  a  set  of  diffuse  functions  built  up  by  first  placing  a  diffuse  p- function  (a=0.018®)  at  the  bond 
center.  To  this  was  added  a  single  diffuse  s-function  whose  exponent  was  chosen  to  optimize  the 
energy  of  the  l*n„  state.  The  diffuse  s-  and  p-functions  were  then  split  into  two  uncontracted 
functions  using  the  prescription  of  Dunning  (1977),  and  to  this  was  added  a  single-component 
d-function  (Peric,  Buenker,  and  Peyerimhoff  1984),  producing  a  (2s,  2p,  Id)  set  of  diffuse  orbitals. 

Each  set  of  d-fimctions  contain  six  angular  momentum  components,  thus  giving  a  total  of  84  basis 
functions  (details  are  given  in  Table  1). 

The  R(C-H)  bond  lengths  were  held  constant  at  the  experimentally  determined  ground  state  value 
of  1.058X  (Herzberg  1966).  Henceforth  all  quantities  shall  be  reported  in  atomic  units  unless  stated 
otherwise.  The  PESs  were  calculated  for  the  )  and  Q'llu )  electronic  states  at  the  points 
R(C-0=(1.783,  1.983,  2.083,  2.183,  2.283,  2.383,  2.483,  2.583,  2.783)  in  the  reduced  symmetry  of  the 
Dj^  point  group  (both  states  have  linear  conformations)  with  the  C-C  bond  as  the  z-axis.  This  section 
of  the  PES  should  be  sufficient  for  estimating  the  dependence  of  the  electric  transition  dipole  moment 
on  nuclear  motion. 
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Table  1.  Atomic  Basis  Set  (McLean  and  Chandler) 


Carbon 

Exponent  C.C. 

Hydrogen 

Exoonent  C.C. 

Diffuse 

Extwnent  C.C. 

Is 

Is 

Is 

15469.40 

0.000570 

82.63637 

0.002006 

0.0304  1.0 

2316.470 

0.004388 

12.40956 

0.015345 

527.099 

0.022810 

2.823854 

0.075577 

149.438 

0.091563 

0.797670 

0.256922 

48.8562 

0.290083 

17.62090 

0.679184 

2s 

2s 

2s 

6.810820 

0.637540 

0.258053 

1.0 

0.012  1.0 

2.7276 

0.397545 

3s 

3s 

0.756740 

1.0 

0.089891 

1.0 

0.300730 

1.0 

0.114090 

1.0 

Ifi 

l£ 

is 

51.72330 

0.009119 

0.75 

1.0 

0.0342  1.0 

12.33970 

0.063253 

3.772240 

0.269747 

1.32487 

0.759674 

2b 

2b 

0.505460 

1.0 

0.0135  1.0 

3B 

0.198270 

1.0 

0.077310 

1.0 

id 

Id 

0.35 

1.0 

0.015  1.0 

2d 

1.50 

1.0 
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The  SA-CASSCJF  calculations  were  perfonned  at  each  point  on  the  PESs  in  order  to  generate  a 
flexible  set  of  oibitals  for  use  in  the  Cl  calculations  (Diffenderfer  and  Yartcony  1981).  The  Cl 
calculations  were  perfonned  with  the  ALCHEMY  codes  of  Liu  and  Yoshimine  (1981).  A  set  of 
reference  CSFs  were  used  to  generate  the  final  Cl  expansions  which  consist  of  all  single  and  double 
excitations  relative  to  these  reference  CSFs.  Further  descriptions  of  the  SA-CASSCF  and  Cl 
calculations  will  be  presented  in  the  appropriate  sections. 

3.1  State  Description.  Since  only  singlet  spin  states  are  discussed  in  this  study,  the  spin  labels 
will  be  dropped  from  the  state  descriptions.  In  the  ground  state,  the  MO  occupations  are 
laglo^^o^o^OglJC^jt^y.  The  components  of  the  C  state  are  represented  by  the  MO  fillings 
...ijt^^Tt^Og  and  ,  where  the  40g  represents  a  s-Rydberg  orbital.  Based  on  these  state 

descriptions,  the  active  orbitals  chosen  for  use  in  the  SA-CASSCF  consisted  of  the  Ibj^), 
ltCuyGb2Q).40g(4ag),  with  the  4iCy  electrons  distributed  in  all  possible  ways  amongst  these  orbitals.  The 
weights,  Wg,  in  Equation  1  were  chosen  to  be  (2,1,1)  for  X,  QIlu,),  and  Ql^y),  in  that  order.  The 
resulting  orbitals  were  then  used  as  expansion  vectors  in  the  Cl  calculations. 

The  reference  CSFs  used  for  generating  the  Cl  wavefimctions  are  listed  in  Table  2.  The  sum  of 
the  squares  of  the  Cl  coefficients,  i.e.,  Zt^j,  for  the  reference  set  varies  from  0.89  to  0.92  for  X  and 
0.90  to  0.93  for  each  C  component  over  the  range  of  the  PESs.  Both  the  and  n„y  components 
were  calculated  as  an  internal  check  on  the  symmetry  of  the  wavefimctions.  As  an  additional  check, 
the  magnetic  transition  dipole  moment,  <Hj>,  coupling  and  n„y  was  calculated  to  be  0.5,  which  is 
correct  for  a  U  state,  t  The  total  number  of  CSFs  used  to  describe  the  electronic  states  are  29,509 


tThe  actual  expression  representing  the  Ae/E  calculated  in  this  study  is  given  by 

which  differs  from  Equation  5  only  in  the  units  conversion  due  to  the  substitution  of  the  orbital  angular  momentum  operator 
for  the  magnetic  dipole  moment  operator  and  carrying  through  the  multiplication  by  /is,  the  Bohr  magneton.  The  definition 
of~^]  then  becomes 


where  use  has  been  made  of  the  relationship  \i  K1/2)L . 
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Table  2.  Reference  CSFs  Used  in 


IRREP: 

“g 

b.g 

b,u 

b2g 

bag 

bsu 

MO; 

8u 

5g 

Ou 

^gx 

^uy 

^gy 

IRREP 

2o„ 

3o„ 

27tu, 

2 

2 

2 

- 

2 

- 

2 

2 

2 

- 

2 

- 

2 

2 

2 

- 

- 

- 

2 

2 

2 

- 

2 

- 

2 

2 

2 

- 

2 

- 

2 

2 

2 

- 

1 

- 

2 

2 

2 

- 

2 

- 

2 

2 

2 

- 

1 

- 

2 

2 

2 

- 

2 

- 

B,u(nj 

2 

2 

1 

2 

- 

1 

- 

2 

2 

1 

2 

- 

1 

- 

BMy) 

2 

2 

1 

2 

- 

2 

- 

2 

2 

1 

2 

- 

- 

- 

B2.(n^) 

2 

2 

2 

1 

1 

- 

2 

2 

2 

1 

I 

- 

B,g(n„) 

2 

2 

2 

1 

2 

- 

2 

2 

2 

1 

- 

- 

Bu(Ag) 

2 

2 

2 

- 

1 

- 

2 

2 

2 

- 

2 

1 

2 

2 

2 

- 

1 

- 

2 

2 

2 

- 

2 

- 

471.,, 


1 


l7C„y  2k^  4k^  iTlp  271, 

2  -  -  - 

2  -  -  2 

1  -  -  1 

1  -  -  -  1 

2  -  -  1 

1  .  .  - 

2  -  -  -  1 

1  .  .  - 

2  -  -  - 


1  -  - 

1  -  -  2 

2  -  - 


1 

1 

2  1 
1 
2 
1 


2 


iTCp,  27Cjy 


1 


1 


2 


2 


1 


for  the  X,  and  96,260  for  each  component  of  the  C  state.  The  PESs  resulting  ftom  the  CIs  are  shown 
in  Figure  1  and  the  values  given  explicitly  in  Table  3. 

3.2  Vibrational  Treatment  Effects  of  nuclear  motion  were  included  using  an  approximate 
symmetric  stretch  normal  mode.  For  all  the  PESs  the  C-H  bond  lengths  were  held  constant,  and  for 
the  C-C  stretch  the  vibrational  wavefunctions  were  calculated  by  treating  the  molecule  as  a 
homonuclear  diatomic.  Two  different  estimates  were  used  for  the  effective  reduced  mass  for  this 
mode.  In  the  first  approach,  atomic  masses  of  13  were  centered  at  each  carbon  atom,  giving  a  reduced 
mass  of  6.5  amu  for  both  states.  In  the  second  approach,  for  each  state  the  effective  reduced  mass  of 
the  system  was  taken  to  reproduce  the  experimental  stretching  fundamental,  Vj.  This  predicts  reduced 
masses  of  X(7.15  amu)  and  Q6.92  amu).  The  differences  in  reduced  masses  result  in  changes  in  the 
final  MCD  intensities  from  2%  in  the  0-0  transition  to  5%  in  the  0-3  transition  which  are  certainly 
acceptable  in  the  context  of  these  calculations.  All  the  vibrationally-averaged  quantities  reported  in 
this  T04850  study  were  obtained  using  the  latter  approach  to  approximating  the  nuclear  masses. 

4.  RESULTS 

The  results  from  the  vibrational  analyses  on  the  X  and  C  states  are  summarized  in  Table  5, 
including  the  energies  for  the  lowest  four  vibrational  states.  Gedanken  and  Schnepp  assign  the 
progression  beginning  at  65,860  cm'*  in  their  MCD  and  absorption  spectra  to  the  totally  symmetric 
C-C  stretch.  This  is  to  be  compared  to  the  calculated  T^  of  63,175  cm'*. 


Figure  2  includes  the  electric  transition  dipole  moment  as  a  function  of  R(C-C)  and  Table  4  lists 
these  values.  These  are  the  first  high  quality  calculations  reported  for  these  quantities  as  a  function  of 
C-C  distance.  The  oscillator  strengths  obtained  from  vibrationally  averaging  these  moments  are  given 
in  Table  6.  The  sum  of  the  f  values  for  the  first  four  vibrational  bands  is  0.157.  An  early 
experimental  detennination  of  the  oscillator  strength  by  Moe  and  Duncan  (1952)  gave 
f  ^*P=0.0616±.(X)92.  They  qualify  their  number  by  stating  that  the  experimental  values  might  best  be 
treated  as  a  minimum  value,  while  "...the  true  values  may  be  much  larger."  A  more  recent 
experimental  value  of  f  ®*p=0.110  was  extracted  from  the  absorption  spectrum  of  Gedanken  and 
Schnepp.  With  this  new  value,  the  comparison  between  theory  and  experiment  has  improved 
considerably,  but  a  discrepancy  of  43%  still  remains  between  theory  and  experiment. 
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Energy  |au| 


R(C-C)  laul 


Figure  1.  Sections  of  the  Potential  Energy  Surfaces  as  a  Function  of  C-C  Distance  for  the 
and  Cd'H.)  States  (All  Values  in  au).  The  C-H  Bond  Lengths  are  Held 
Fixed  at  the  Ground  State  Value. 
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Table  3.  X  and  C('nj  State  Energies  as  a  Function  of  C-C  Bond  Length 


R(C-0 

au 

X(‘Zg*) 

au 

aU) 

au 

1.78279 

-76.91684 

-76.56676 

1.98279 

-77.07607 

-76.75395 

2.08279 

-77.11386 

-76.80332 

2.18279 

-77.13276 

-76.83253 

2.28279 

-77.13758 

-76.84657 

2.38279 

-77.13202 

-76.84926 

2.48279 

-77.11889 

-76.84356 

2.58279 

-77.10039 

-76.83175 

2.78279 

-77.05356 

-76.79637 

Table  4.  Electric  Transition  Dipole  Moment  as  a  Function  of  R(C-C) 


R(C-0 

au 

au 

1.7828 

0.4983 

1.9828 

0.5616 

2.0828 

0.5891 

2.1828 

0.6139 

2.2828 

0.6362 

2.3828 

0.6560 

2.4828 

0.6734 

2.5828 

0.6886 

2.6828 

0.7129 

♦  ***  '^ 

The  magnitude  of  the  X  component  is  identical  to  at  least  the  number  of  digits  reported. 
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Table  5.  Results  of  the  Vibrational  Analysis  on  the  X  and  C  States  for  the 
Symmetric  C-C  Stretch 


Predicted  Vibrational  Energies  fem'^) 

V 

X 

C 

0 

988 

921 

1 

2955 

2749 

2 

4905 

4563 

3 

6829 

6369 

■Uj  Fundamental* 

X 

C 

Theory 

T 

*0 

63175 

(7.83eV) 

Theory 

1977 

1841 

Exper. 

65814 

(8.16eV)'’ 

Exper.** 

1974 

1840 

Exper. 

65860 

(8.17eV)‘= 

a.  The  theoretical  U2  obtained  from  varying  the  reduced  mass  of  the  system,  and  is  defined  as 
‘U2=2  X  E(v=0). 

b.  From  Herzbcrt  (1966). 

c.  From  Gedanken  and  Schnepp  (1976). 


Table  6.  MCD  Intensities  and  Oscillator  Strengths  for  Vibrational  Bands  in  the  C<“X  Transition 


v" 

D 

<A£>,^ 

f  Cdc 

f 

f  BV* 

f  TlwMy® 

0 

0 

0.0179 

0.0125±.001 

0.1046 

0.0743 

0 

1 

0.0068 

0.00401.0004 

0.0408 

0.0306 

0 

2 

0.0015 

0.00071.0001 

0.0094 

0.0046 

0 

3 

0.0003 

0.00011.00001 

0.0017 

0.0008 

TOTAL 

0.0265 

0.0173 

0.1565 

0.1103 

0.0616  .0092 

0.072 

a)  Reduced  mass  determined  by  fitting  calculated  V2  to  experiment  See  text  for  details. 

b)  From  Gedanken  and  Schnepp  (1976). 

c)  Calculated  from  the  2^  values  listed  in  Table  2  of  Gedanken  and  Schnepp  (1976).  These  were  extracted  by  Gedanken 
and  Schnepp  from  their  absorption  spectrum  after  subtracting  out  intensity  from  an  unassigned,  imderlying  band. 

d)  See  Moe  and  Ehmcan  (1952). 

e)  From  ab  initio  "frozen  core"  calculations.  Includes  no  correlation  corrections.  See  Demoulin  and  Jungen  (1974). 


13 


The  MCD  spectrum  provides  another  independent  check  on  the  electric  transition  dipoles.  The 
calculated  and  experimental  MCD  intensities  for  the  first  four  vibrational  bands  are  also  given  in 
Table  6.  The  agreement  is  qualitatively  correct,  showing  decreasing  MCD  intensity  with  the 
increasing  v'.  One  finds  reasonably  good  agreement  between  experiment  and  theory  for  the  0-0 
transition  where  theory  and  experiments  predict  <A£>  ,=0.0179  and  0.0125±.001,  respectively.  Again, 
the  theory  predicts  an  intensity  that  is  43%  larger  than  experiment;  consistent  with  the  difference  seen 
in  absorption.  The  reader  should  be  reminded  that  the  same  electric  dipole  moments  are  used  in  the 
theoretical  predictions  of  the  oscillator  strength  and  the  MCD  intensities.  This  consistency  in  the 
comparisons  between  theory  and  experiment  for  absorption  intensities  and  MCD  intensities  supports 
the  experimental  correction  to  the  absorption  intensity  achieved  by  subtracting  out  the  broad  band 
underlying  the  experimental  absorption  spectrum.  The  values  for  the  1-0,  2-0,  and  3-0  bands  are 
(theory,experiment)=(0.0068,0.0040),o,  (0.0015,0.0007)2o.  (0.0003 ,0.000 1)30.  All  the  calculated  values 
are  seen  to  be  larger  than  their  experimental  counterparts.  Assuming  the  experimental  MCD  values  are 
accurate,  the  fact  that  the  lower  v'  levels  agree  better  with  experiment  than  the  higher  v'  levels 
suggests  increasing  inaccuracies  in  the  Franck-Condon  overlap  with  these  higher  vibrational  levels. 

This  in  turn  would  indicate  that  the  C  state  PES  is  not  described  with  equal  accuracy  over  the  range  of 
the  PES  section  studied.  Specifically,  the  PES  seems  better  described  near  the  equilibrium,  R(C-C). 

5.  CONCLUSIONS 

This  is  the  first  report  using  a  large  atomic  basis  set  and  large  Cl  expansions  to  predict  a  priori 
the  electric  dipole  transition  moment  for  the  C^X  band  system  as  a  function  of  R(C-C).  The 
calculated  To=63,175  cm'*  is  in  reasonable  agreement  with  the  experimental  value  of  65,860  cm'*. 

An  approximate  vibrational  averaging  of  the  electric  transition  dipole  moment  predicts  an  absorption 
oscillator  strength  of  f=0.157,  which  is  more  than  a  factor  of  2  larger  than  an  early  experimental  value 
of  f=0.062±0092,  but  in  better  agreement  with  a  more  recent  value  of  0.1 10  extracted  from  the 
experimental  absorption  spectrum  of  Gedanken  and  Schnepp  (1976). 

The  experimental  MCD  intensities  for  the  0“X  transition  was  also  used  as  another  check  on  the 
calculated  electric  transition  dipole  moment.  The  advantage  of  using  the  MCD  over  the  absorption  is 
the  apparent  lack  of  interference  in  the  band  structure  from  a  broad  band  present  in  the  absorption. 

The  integrated  MCD  intensity,  <Ae>,,  for  the  0-0  transition  is  predicted  to  be  0.0179,  which  can  be 
compared  to  the  experimental  value  of  0.0125±.(X)1.  Theory  also  predicts  the  correct  trend  of 


14 


decreasing  MCD  intensity  for  transitions  to  the  higher  v'  levels  of  C  even  though  the  difference 
between  theory  and  experiment  for  <Ae>|  increases  with  increasing  v'.  The  remaining  discrepancy  of 
approximately  40%  between  theory  and  experiment  in  both  the  MCD  and  absoiption  intensities 
suggests  the  need  for  more  experimental  work. 
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2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  R.  Bemecker,  R-13 

G.B.  Wilmot,  R-16 
Silver  Spring,  MD  20903-5000 
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5  Commander 

Naval  Research  Laboratory 
ATTN:  M.C.  Lin 

J.  McDonald 
E.  Oran 
J.  Shnur 

RJ.  Doyle,  Code  6110 
Washington,  DC  20375 

1  Commanding  Officer 
Naval  Underwater  Systems 

Center  Weapons  Dept 
ATTN:  R.S.  Lazar/Code  36301 

Newport,  RI  02840 

2  Commander 

Naval  Weapons  Center 
ATTN:  T.  Boggs,  Code  388 

T.  Parr,  Code  3895 
China  Lake,  CA  93555-6001 

1  Superintendent 

Naval  Postgraduate  School 
Dept,  of  Aeronautics 
ATTN:  D.W.  Netzer 

Monterey,  CA  93940 

3  AULSCF 

ATTN:  R.  Corley 

R.  Geisler 
J,  Levine 

Edwards  AFB,  CA  93523-5000 

1  AL/MKPB 

ATTN:  B.  Goshgarian 

Edwards  AFB,  CA  93523-5000 

1  AFOSR 

ATTN:  J.M.  Tishkoff 

Bolling  Air  Force  Base 
Washington,  DC  20332 

1  OSD/SDIO/IST 
ATTN:  L.  Caveny 

Pentagon 

Washington,  DC  20301-7100 

1  Commandant 
USAFAS 

ATTN:  ATSF-TSM-CN 

Fort  Sill,  OK  73503-5600 

1  FJ.  Seiler 

ATTN:  Sj\.  Shackleford 

USAF  Academy,  CO  80840-6528 
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1  University  of  Dayton  Research  Institute 
ATTN:  D.  Campbell 

AIVPAP 

Edwards  AFB,  CA  93523 
1  NASA 

Langley  Research  Center 
Langley  Station 

ATTN:  G.B.  Northam/MS  168 

Hampton,  VA  23365 

4  National  Bureau  of  Standards 
ATTN:  J.  Hastie 

M.  Jacox 
T.  Kashiwagi 
H.  Semerjian 

US  Department  of  Commerce 
Washington,  DC  20234 

1  Aerojet  Solid  Propulsion  Co. 

ATTN:  P.  Micheli 

Sacramento,  GA  95813 

1  Applied  Combustion  Technology,  Inc. 

ATTN:  A.M.  Varney 

P.O.  Box  607885 
Orlando,  FL  32860 

2  Applied  Mechanics  Reviews 
The  American  Society  of 

Mechanical  Engineers 
ATTN:  R.E.  White 

A.B.  Wenzel 
345  R  47th  Street 
New  York,  NY  10017 

1  Atlantic  Research  Corp. 

ATTN:  M.K.  King 

5390  Cherokee  Avenue 
Alexandria,  VA  22314 

1  Atlantic  Research  Corp. 

ATTN:  R.H.W.  Waesche 

7511  Wellington  Road 
Gainesville,  VA  22065 

1  AVCO  Everett  Research 
Laboratory  Division 
ATTN:  D,  Stickler 

2385  Revere  Beach  Parkway 
Everett,  MA  02149 


No.  of 

Copies  Organization 

1  Battelle  Memorial  Institute 
Tactical  Technology  Center 
ATTN:  J.  Huggins 

505  King  Avenue 
Columbus,  OH  43201 

1  Cohen  Professional  Services 
ATTN:  N.S.  Cohen 

141  Channing  Street 
Redlands,  CA  92373 

1  Exxon  Research  &  Eng.  Co. 
ATTN:  A.  Dean 

Route  22E 

Annandale,  NJ  08801 

1  Ford  Aerospace  and 
Communications  Corp. 
DIVAD  Division 
Div.  Hq.,  Irvine 
ATTN:  D.  Williams 

Main  Street  &  Ford  Road 
Newport  Beach,  CA  92663 

1  General  Applied  Science 
Laboratories,  Inc. 

77  Raynor  Avenue 
Ronkonkama,  NY  11779-6649 

1  General  Electric  Ordnance 

Systems 

ATTN:  J.  Mandzy 

100  Plastics  Avenue 
Pittsfield,  MA  01203 

2  General  Motors  Rsch  Labs 

Physics  Department 
ATTN:  T.  Sloan 

R.  Teets 
Warren.  MI  48090 

2  Hercules,  Inc. 

Allegheny  Ballistics  Lab. 

ATTN:  W.B.  Walkup 

E.A.  Yount 
P.O.  Box  210 

Rocket  Center,  WV  26726 

1  Honeywell,  Inc. 

Government  and  Aerospace 
Products 

ATTN:  D.E.  Broden/ 

MS  MN50-2000 
600  2nd  Street  NE 
Hopkins,  MN  55343 
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1  Honeywell,  Inc. 

ATTN:  R.E.  Tompkins 

MN38-3300 

10400  Yellow  Circle  Drive 
Minnetonka,  MN  55343 

1  IBM  Corporation 
ATTN:  A.C.  Tam 

Research  Division 
5600  Cottle  Road 
San  Jose,  CA  95193 

1  in  Research  Institute 

ATTN:  R.F.  Remaly 

10  West  35th  Street 
Chicago,  IL  60616 

2  Director 
Lawrence  Livermore 

National  Laboratory 
ATTN:  C  Westbrook 

M.  Costantino 
P.O.  Box  808 
Livermore,  CA  94550 

1  Lockheed  Missiles  &  Space  Co. 
ATTN;  George  Lo 
3251  Hanover  Street 
Dept.  52-35/B204/2 
Palo  Alto,  CA  94304 

1  Los  Alamos  National  Lab 
ATTN:  B.  Nichols 

T7,  MS.B284 
P.O.  Boot  1663 
Los  Alamos,  NM  87545 

1  National  Science  Foundation 
ATTN:  A.B.  Harvey 

Washington,  DC  20550 

1  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 

P.O.  Box  222 

St.  Marks,  FL  32355-0222 

1  Paul  Gough  Associates,  Inc. 
ATTN:  P.S.  Gough 

1048  South  Street 
Portsmouth,  NH  03801-5423 


No.  of 
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2  Princeton  Combustion 

Research  Laboratories,  Inc. 
ATTN:  M.  Summerfield 

NA.  Messina 
475  US  Highway  One 
Monmouth  Junction,  NJ  08852 

1  Hughes  Aircraft  Company 
ATTN:  T.E.  Ward 

8433  Fallbrook  Avenue 
Canoga  Park,  CA  91303 

1  Rockwell  International  Corp. 
Rocketdyne  Division 
ATTN:  J.E.  F1anagan/HB02 

6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

4  Sandia  National  Laboratories 
Division  8354 
ATTN:  R,  Cattolica 

S.  Johnston 
P.  Mattem 
D.  Stephenson 
Livermore,  CA  94550 

1  Science  Applications,  Inc. 

ATTN:  R.B.  Edelman 

23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

3  SRI  International 
ATTN;  G.  Smith 

D.  Crosley 
D.  Golden 

333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of  Tech. 
Davidson  Laboratory 
ATTN:  R.  McAlevy,  III 

Hoboken.  NJ  07030 

1  Sverdrup  Technology,  Inc. 

LeRC  Group 

ATTN:  RJ.  Locke,  MS  SVR2 

2001  Aerospace  Parkway 
Brook  Park,  OH  44142 

1  Thiokol  Corporation 
Elkton  Division 
ATTN:  S.F.  Palopoli 

P.O.  Box  241 
Elkton,  MD  21921 
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1  Morton  Thiokol,  Inc. 

Huntsville  Division 
ATTN:  J.  Deur 

Huntsville,  AL  35807-7501 

3  Tbiokol  Corporation 
Wasatch  Division 
ATTN:  SJ.  Bennett 

P.O.  Box  524 
Brigham  City,  UT  84302 

1  United  Technologies  Research  Center 
ATTN:  A.C.  Eckbreth 

East  Hartford,  CT  06108 

3  United  Technologies  Corp. 

Chemical  Systems  Division 
ATTN:  R.S.  Brown 

T.D.  Myers  (2  copies) 

P.O.  Box  49028 

San  Jose,  CA  95161-9028 

1  Universal  Propulsion  Company 
ATTN:  HJ.  McSpadden 

Black  Canyon  Stage  1 
Box  1140 

Phoenix,  AZ  85029 

1  Veritay  Technology,  Inc. 

ATTN:  RB.  Fisher 

4845  Millersport  Highway 
P.O.  Bok  305 

East  Amherst,  NY  14051-0305 


1  University  of  California, 
Berkeley 

Chemistry  Deparment 
ATTN:  C.  Bradley  Moore 
21 1  Lewis  Hall 
Berkeley,  CA  94720 

1  University  of  California, 

San  Diego 

ATTN:  FA-  Williams 

AMES,  BOlO 
La  Jolla,  CA  92093 

2  University  of  California, 

Santa  Barbara 
Quantum  Institute 
ATTN:  K.  Schofield 

M.  Steinberg 

Santa  Barbara,  CA  93106 

1  University  of  Colorado  at 

Boulder 

Engineering  Center 
ATTN:  J.  Daily 

Campus  Box  427 
Boulder,  CO  803094M27 

2  University  of  Southern 

California 
Dept,  of  Chemistry 
AT^:  S.  Benson 

C.  Wittig 

Los  Angeles,  CA  90007 


1  Brigham  Young  University 
Dept,  of  Chemical  Engineering 
ATTN:  M.W.  Beckstead 

Provo,  UT  84058 

1  California  Institute  of  Tech. 

Jet  Propulsion  Laboratory 
ATTN:  L.  Strand/MS  512/102 

4800  Oak  Grove  Drive 
Pasadena,  CA  91009 


1  Case  Western  Reserve  Univ. 
Div.  of  Aerospace  Sciences 
ATTN;  J.  Tien 
Cleveland,  OH  44135 

1  Cornell  University 

Department  of  Chemistry 
ATTN:  TA.  Cool 

Baker  Laboratory 
Ithaca,  NY  14853 


1  California  Institute  of 
Technology 

ATTN:  F.E.C.  Culick/ 

MC  301-46 
204  Karman  Lab. 

Pasadena,  CA  91125 

1  University  of  California 
Los  Alamos  Scientific  Lab. 

P.O.  Box  1663,  Mail  Stop  B216 
Los  Alamos,  NM  87545 


1  University  of  Delaware 
ATTN:  T.  Brill 

Chemistry  Department 
Newark,  DE  19711 

1  University  of  Florida 
Dept,  of  Chemistry 
ATTN:  J.  Winefordner 

Gainesville,  FL  32611 
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3  Georgia  Institute  of 
Technology 
School  of  Aerospace 
Engineering 
ATTN:  E.  Price 

W.C.  Strahle 
B.T.  Zinn 
Atlanta,  GA  30332 

1  University  of  Illinois 
Dept  of  Mech.  Eng. 

ATTN:  H.  Krier 

144MEB,  1206  W.  Green  St 
Urbana,  IL  61801 

1  Johns  Hopkins  University/APL 
Chemical  Propulsion 
Information  Agency 
ATTN:  T.W.  Christian 

Johns  Hopkins  Road 
Laurel,  MD  20707 

1  University  of  Michigan 
Gas  Dynamics  Lab 
Aerospace  Engineering  Bldg. 
ATTN:  G.M.  Facth 

Ann  Arbor,  MI  48109-2140 

1  University  of  Minnesota 
Dept  of  Mechanical 
Engineering 
ATTN:  E.  Fletcher 

Minneapolis,  MN  55455 

3  Pennsylvania  State  University 
Applied  Research  Laboratory 
ATTN:  ICK.  Kuo 

H.  Palmer 
M.  Micci 

University  Park,  PA  16802 

1  Pennsylvania  State  University 
Dept,  of  Mechanical  Engineering 
ATTN:  V.  Yang 

University  Park,  PA  16802 

1  Polytechnic  Institute  of  NY 
Graduate  Center 
ATTN:  S.  Lederman 

Route  110 

Farmingdale,  NY  11735 


2  Princeton  University 
Forrcstal  Campus  Library 
ATTN:  K.  Brezinsky 

L  Glassman 
P.O.  Box  710 
Princeton,  NJ  08540 

1  Purdue  University 
School  of  Aeronautics 
and  Astronautics 
ATTN:  J.R.  Osborn 

Grissom  Hall 

West  Lafayette,  IN  47906 

1  Purdue  University 

Department  of  Chemistry 
ATTN:  E.  Grant 

West  Lal^yette,  IN  47906 

2  Purdue  University 
School  of  Mechanical 

Engineering 

ATTN:  N.M.  Laurendeau 

S.N.B.  Murthy 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic  Inst. 
Dept,  of  Chemical  Engineering 
ATTN:  A.  Fontijn 

Troy,  NY  12181 

1  Stanford  University 
Dept,  of  Mechanical 
Engineering 
ATTN:  R.  Hanson 

Stanford,  CA  94305 

1  University  of  Texas 
Dept,  of  Chemistry 
ATTN:  W.  Gardiner 

Austin,  TX  78712 

1  University  of  Utah 

Dept,  of  Chemical  Engineering 
ATTO:  G.  Flandro 

Salt  Lake  City,  UT  84112 

1  Virginia  Polytechnic 
Institute  and 
State  University 
ATTN:  J.A.  Schetz 

Blacksburg,  VA  24061 
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1  Freedman  Associates 
ATTN:  E.  Freedman 

2411  Diana  Road 
Baltimore,  MD  21209-1525 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratofy  undertakes  a  continuing  effort  to  improve  the  ^ality  of  the  reports  it  publishes. 
Your  comments/answers  to  die  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  brl-tr-3183 _ Date  of  Report  December  1990 

2.  Date  Report  Received _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _  .  _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so.  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  formtu,  etc.)  . . 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  Suite,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21(>i  -5066 
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t 
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IN  THE  i 

UHRED  STATES  | 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PER^y!lT  No  00C1,  APG,  MD 


POSTAGE  Will  P£  PAID  BY  AODRcSSEc 


Director 

U.S.  A.nny  BulUsiic  Rcsu'ah  Lah-oralory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MD  21005-9589 


FOLD  HERE 


